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Abstract 

For two-black-hole time-symmetric initial data we consider the Corvino con- 
struction of gluing an exact Schwarzschild end. We propose to do this by using 
Brill waves. We address the question of whether this method can be used to re- 
duce the overall energy, which seems to relate to the question of whether it can 
reduce the amount of 'spurious' gravitational radiation. We find a positve answer 
at first order in the inverse gluing radius. 

1 Introduction 

Processes involving collisions and mergers of black holes are amongst the most 
promising candidates regarding the generation of sufficiently strong gravitational 
waves. Hence much analytical and numerical effort went into the modeling of such 
processes as Cauchy problems for the matter-free Einstein equations. Initial data 
consist of a triple (a, fi, K), where L is a three manifold (the Cauchy hypersurface) 
and h, K are two symmetric second-rank tensor fields satisfying a system of four 
underdetermined-elliptic differential equations, the constraints, h is a Riemannian 
metric for L and K essentially its time derivative. Here we shall only be concerned 
with data representing two black holes at the moment of time symmetry. The latter 
condition is equivalent to K = 0. Physically this corresponds to the most simple sit- 
uation where two uncharged and unspinning black holes are placed at some distance, 
without relative velocity and without any overall linear or angular momentum. We will 
even assume equal hole masses so that the only free parameters are the mass and the 
mutual distance. 

Data representing two holes at the moment of time symmetry have been known 
for long and can be written down analytically (see e.g. |5 ] for a recent review). One 
might think that the really hard part is their time evolution by means of a sophisticated 
combination of analytical and numerical techniques. However, there are also unsolved 
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conceptual and analytical problems regarding just the initial data. One of them is the 
problem of 'spurious' gravitational radiation, which we address here. 



V 



4) 




I 



Figure 1: Conformal diagram of spacetime with two spacelike hypersurfaces L and 
L' ending at spacelike infinity io and future lightlike infinity J + respectively. The 
difference between the ADM-energy of L (computed at io) and the Bondi mass of 1 ' 
(computed at the intersection 2-sphere between L' and 3 + , here denoted by the point 
b) must be due to gravitational radiation escaping between io and b. This radiation 
originates form the causal past P of the region iob C J + , whose intersection with L 
is the region aio- In this sense we say that the data on L contained radiation in that 
region. Any gravitational radiation emerging from an (quasi localized) event e, e.g. 
the formation of a black hole due to neutron-star mergers, cannot register before c, the 
intersection of the future light-cone at e with 3 + . 

It has recently been shown to some level of rigour that even the most simple two 
black-hole collision data contain radiation in any neighbourhood of spacelike infinity 
i°; see Figure^ This was done by approximately evaluating the Bondi mass on a cut 
of J + close to i° for time symmetric, conformally flat data. The result was [ 8 ] 



where the are the Newman-Penrose constants [6] and where u is the Bondi pa- 
rameter along J + , which tends to — oo as one approaches i°. Hence m-Bondi < ht-adm 
with equality at i , as it must be. If some of the Newman-Penrose constants do not 
vanish, we have a strict inequality off i°, meaning that some gravitational radiation 
has reached 3 + in between. If this is the case, one says that the initial data contained 
radiation in a neighbourhood of (spacelike) infinity. This radiation is spurious in the 
sense that it cannot be caused by the collision process. Hence one would like to re- 
move it. If we were in a linear theory, e.g. electrodynamics, the meaning of 'remove' 
would be to subtract the appropriate solution to the homogeneous equations. However, 
in a non-linear theory, like General Relativity, it is a priori unclear how to analytically 
achieve a 'removal'. Is it possible at all to get rid of the spurious radiation without 
touching the interior data? 



ITLBondi = ttladm + Y |G k | 2 — + 0(l/u 8 ) 
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The question addressed here relates to the problem of analytically understanding 
the region where null infinity touches spacelike infinity. The first step in this direction 
was taken by Friedrich [4|, who used the conformal field equations to analyse this re- 
gion. In this approach, spacelike infinity is blown up to a cylinder S 2 x [0, 1], whose two 
boundary components connect to future and past null infinity (J ± ) respectively. The 
Cauchy surface, L, intersects the cylinder in a 2-sphere cross section, i , from which 
certain transport equations govern the propagation of initial data along the cylinder 
to its boundary. This propagation will generically lead to logarithmic divergences at 
these intersections and, consequently, to a non-smooth null infinity. Whether, and to 
what degree, a non-smooth null infinity is physically objectionable is currently unclear. 
Some necessary conditions for the avoidance of such singularities are known J9j, but 
so far no theorem establishes a sufficient set. Instead, the following conjecture was 
launched [9|: 

Conjecture: If an initial data set which is time-symmetric and confomally flat in a 
neighbourhood B a (io) of the point io yields a development with a smooth null infinity, 
then the initial data are, in fact, Schwarzschild in B Q (i). 

If true, the conjecture tells us that gravitational radiation is excluded near space- 
like infinity. At this point one might fear that the condition of being asymptotically 
Schwarzschild might exclude most data of physical interest. But fortunately this is 
not true, thanks to the work of Corvino [31, who proved that is possible to glue a 
Schwarzschild metric along an annulus to an asymptotically flat, conformally flat 
three metric if the gluing radius is chosen sufficiently large and if the mass of the 
Schwarzschild metric is chosen appropriately. That is, there do exist spacetimes which 
are Schwarzschild at infinity (therefore ensuring a smooth null-infinity) and non-static 
in the inside. 

For example, we could glue a Brill-Lindquist initial data set for the collision of 
two black holes to a Schwarzschild metric. This would clearly remove the 'spurious' 
radiation at spacelike infinity. However, it might at the same time introduce 'spurious' 
radiation through the modified gluing metric on the annulus, which is generally very 
complicated. In particular, it is not conformally flat. This motivates to ask about the 
mass of the glued Schwarzschild solution, depending on the specific gluing instruc- 
tions. 

Information about the mass is not easy to obtain: the gluing function is hard to 
control and not easily calculated explicitly. What we need is a physical idea of how 
the gluing function might look like. Here we propose to specialize to Brill waves. 
They are axisymmetric (as is Brill-Lindquist) and not conformally flat. Is it possible 
to construct a spacetime which is Brill-Lindquist to some radius, has Brill-form on an 
annulus, and is Schwarzschild outside? It is the purpose of this paper to explore this 
idea a little further, based on the diploma thesis by one of us 171. 

2 Corvino 's construction 

We mentioned that the Brill-Lindquist data seem to contain spurious gravitational ra- 
diation near spacelike infinity and, intimately related to this fact, cannot lead to an 
asymptotically simple spacetime. Hence we wish to modify them near infinity. The 
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techniques of Corvino use the underdetermined ellipticity of the constraint equations 
to glue a Schwarzschild end to an arbitrary asymptotically flat manifold. 
Consider a smooth gluing function (3 > 0, | (3 1 < 1 , where 

f 1 forx<1 
P \ forx>2 {Z) 

and an arbitrary asymptotically flat metric gAF together with a Schwarzschild metric 
9ss' C w * tn ADM-energy m and center c. Then, in an asymptotically flat chart, the 
metric 

99 lue = P (^) 9 AF + (l " P (^)) 9F/ 0) 

(for constant R) is smooth and glues g af to g^ c - However, the constraint equa- 
tion of vanishing scalar curvature, R(g g i U e) = 0, is now violated on the annulus, 
R < |x| < 2R, hence g g vue is not a time-symmetric initial datum for Einstein's equa- 
tion. Corvino now shows that it is possible to choose the gluing radius R large enough 
to find a metric H, having support only on the annulus, such that R(g g i U e + H) = 
everywhere. To make this construction work, the variables m and c have to be chosen 
appropriately. In this way Corvino's construction provides us with time-symmetric 
initial data, isometric to Schwarzschild outside some finite radius and non-trivial in 
the inside. But his construction is not very explicit. We have no method at hand to 
construct the metric explicitly, because we do not know much more about the pertur- 
bation h on the annulus than its existence. This leaves important physical questions 
unanswered: Can we give a physical interpretation for the gluing construction? One 
is tempted to ask the following 'physical' questions: does the construction actually 
remove spurious gravitational-radiation energy, or does it at best just lead to a redistri- 
bution of energy from infinity to the annulus region ? The answer does not seem to be 
known. 



3 Brill waves and the gluing construction 
3.1 Brill-Lindquist data 

We propose to phrase these questions in a more explicit environment. Consider Brill- 
Lindquist data describing the head-on collision of two black holes: 

We further simplify by choosing ai = 012, corresponding to equal masses for the holes 
and C] = — C2 = c := (0, 0, d/2), so that the holes lie symmetrically about the origin 
on the z-axis: 

9 « = ( ,+ i?Vi?Tl) 4i> «- <5) 

These data are axisymmetric. We shall equip them with a Schwarzschild end at infinity. 
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3.2 Brill waves 



Our key assumption is that the metric on the annulus is a conformally transformed 
Brill wave fl]. In cylindrical coordinates it correspond to the regular 3-metric 

QBrtii = e 2q(p ' z) (dp 2 + dz 2 ) + p 2 d0 2 . (6) 

Here q is a function of p and z satisfying the following conditions. 

q = if p = (on the z-axis) , 
3 

— q = if p = (on the z-axis) , (7) 
3p 

q — > faster than | for r — > oo . 

The second condition is immediately obvious: If q is to represent an axisymmetric 
solution, it must have an extremum on the axis in the p direction. The first condition 
ensures that there is no conical singularity on the z-axis. The third condition will be 
satisfied automatically because our q will have compact support. 

To solve the constraint equation R = 0, we need to perform a conformal trans- 
formation of the Brill metric. We ask if there exists a conformal factor, ip, with the 
following properties 

R( VWu) = , 

4> > everywhere , (8) 
4> — > 1 at oo . 

Using the behaviour of the Ricci scalar under conformal transformations, R(\b 4 g) = 
— 84>~ 5 (A g — ^R(g)t)j), we end up with the equation 

8A g ib-R(g)il) = 0, (9) 

where \\> > everywhere and \\) — > 1 at oo. In the past much work has been done 
to derive existence theorems for solutions tJj given q. Here we will take a slightly 
different viewpoint. 



3.3 Our construction 

Combining the spirit of Corvino's construction with our physical idea of a Brill wave 
on the annulus, we define a conformal factor 

» + » WWeVfi-pfNeYlf, + (io) 



|x + c|y ' V * / V V R // V I* 

We have generalized the gluing function (3 to be axisymmetric on the annulus: 

f 1 for x < 1 and all 9 
' \ forx > 2 and all ( } 

and (3 (n) (1) = |3 (n) (2) =0foralln > 1. 
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Next we use the fourth power of the function \\) in dTOb as conformal factor to 
define the physical initial-data metric: 

9 = V 9Briii = V (e 2q(p ' z) (dp 2 + dz 2 ) + p 2 d0 2 ) . (12) 

Here q is, by definition, a function with support on the annulus, vanishing deriva- 
tives on its boundary 1 and which satisfies the Brill-conditions @. Taking all this into 
account we observe that 

• g is Brill-Lindquist for r < R, 

• g is Schwarzschild for r > 2R (with ADM-energy M = y), 

• g is a complicated, conformally transformed Brill wave on r G [R, 2R]. 
Now, the constraint equation R(g) = is satisfied if 

( 3 2 3 2 \ , , / 3 2 1 3 1 3 2 \ , ^ -A< 3 'i1j 

( W + ^) q(P ' Z) = \M + r Tr + V- W) q(T > 9) = ' (13) 

Since 4> stands for the expression (flOl . this equation should be read as relating the 
gluing function |3 to the Brill-function q . It turns into an identity (0 = 0) for r < R 
and r > 2R. The mathematical structure of the equation is discussed in the next 
section. 



3.4 The mathematical problem 

Equation (fT3l can be understood as an elliptic equation on the two-dimensional plane. 
However, the boundary conditions are not of the usual Dirichlet or Neumann type. 
Note that the right hand side — for any gluing function |3 — is of bump-type, i.e. van- 
ishes with all derivatives at r = R and r = 2R. Hence we have a large freedom in 
the choice of the right hand side by exploiting the freedom of the gluing function (3, 
of which only the boundary values (including all derivatives) are fixed. On the other 
hand, there are boundary conditions on q and its derivatives that have to be taken into 
account. 

As it turns out, the problem can be formulated in the following way: Given the 
Poisson equation 

A (2) q=f = _^ (14) 

in the half-annular region H sketched in Figure|2 Is it possible to choose the right hand 

side of (fPfl l such that q simultaneously satisfies Dirichlet and Neumann conditions for 

the given region H? This we will call a DN-solution. Note that imposing Dirichlet and 

Neumann conditions guarantees the Brill-conditions (0 on the z-axis, as well as the 

vanishing of all radial derivatives on the arches of H. Indeed, if q and its first derivative 

vanish on the arches, then, by equation (fT4l . all derivatives of q will vanish there since 

f vanishes with all derivatives on the arches. Although one can easily write down some 

f 's admitting a DN-solution, it is not clear if these f can be obtained by — A ^ ^ with 

an appropriate choice of (3 — the root of the problem being that the inhomogeneity is 

highly non-linear in (3. 

'in the sequel a function with these properties will sometimes be called a "function of bump- 
character" (although it can become negative, of course). 
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Figure 2: Half-annular region H (shaded region) in upper z-p halfplane. 



3.5 Approximate results 

Corvino's construction only works for large gluing radii R. We could therefore try to 
solve the DN-problem to first order in the inverse gluing radius. To do this we expand 
the inhomogeneity to first order in 1 /R, making two further choices: 

1 The first is related to the possibility of reducing the energy. The further outside 
the gluing is performed, the smaller the difference of the Schwarzschild and 
the Brill-Lindquist data will be. In other words, the difference will scale in 
some way with the inverse gluing radius — the exact form being dependent on 
the details of the gluing procedure. Recall equation <fT0b . introduce the scaling 
|x| — > |x|R and expand in the inverse gluing radius, ^: 

„ = 1 + Mr , e) (^ + ^ +0 ^)) (15) 

Note that we have no - R ^ z -(dipole)-terms due to the choice of equal masses for the 
two Brill-Lindquist holes and equal mass-centres for the Brill-Lindquist system 
and the Schwarzschild hole. We choose the scaling to be 

2a- A = ^ (16) 

where b takes real values in a small interval around (eventually we will be 
interested in b > 0). 2 

2 The second choice is an ansatz for the gluing function (3 : 

P(t,9) = <x(t) + u(r)sin 2 e. (17) 

Here oc(r) is a radial gluing function (<x(1 ) = 1 , <x(2) = 0, all derivatives vanish 
at these two points) and u(r) is a function with bump character on the annulus 

2 If we chose 2 a — A = ^ we would, at the lowest order, have the same situation as for the gluing 
of two Schwarzschild black holes of ADM-energies 2a and A respectively. On account of the Penrose 
inequality one expects this not to work with mass reduction. Note that for large gluing radii the con- 
formally transformed Brill wave will certainly not introduce a new apparent horizon. In fact, in can be 
shown that the recursion technique presented below fails for the Schwarzschild-Schwarzschild case. 
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= 0, u(2) = 0, all derivatives vanish at these two points) but otherwise 
arbitrary. For the expansion dlOt of \J) we explicitly introduce the gluing radius R 
by scaling r — > Rr in the formulae. This makes the gluing occur on the interval 
[1,2]). 



The right hand side can then be written as 



-^(«'(T) + ?«'iT)+^mr: 



b ad \ , , 

- + ^3 ) + 2a' (r) 



3 ad" 



+ sin 2 



u» + -u'(r) 
r 



b ad 2 
r 2r 3 



3ad 2 
4r 3 



a"fr) + -a'ff 



, , , b 3ad 2 \ 9ad 2 a'(r) 9ad 2 u(r) 
+ 2u (r) ( — T - I + 



2r 4 ; 



2r 4 



+ sin 4 9 



3 ad 2 
4r 3 



r 



9ad 2 u'(r) 6ad 2 



2r 4 



+ 



-u(r) 



(18) 



+ 0(R~ 6 ) 

the right-hand side of which may be written in short-hand: 

^5 (A (a', a", u) + A 2 (a', a", u, u', u") sin 2 
+ A 4 (^Li / ,Li /, )sin 4 0) +0(R 
We make the following ansatz for q(r, 0): 

q(r,0) = -j ^B 1 (r/R)sin 2 + B 2 (r/R)sin 4 

The functions B| are of bump-character in [1,2]. This ansatz satisfies the three Brill- 
conditions © for q(r, 0). The expression will be of order when acted upon by the 
Laplacian. If we apply the Laplacian and perform the scaling with the gluing radius 
afterwards, we obtain the following left hand side: 



(19) 



(20) 



LHS 



1 

R5 



2Bt(t 



+ sin 2 ( B'j(r) + -B\{t) 
+ sin 4 0fB^(r) + lB^r) 



12B 2 (r] 



16 



B 2 W 



r 

The right-hand side is given by minus the expansion ( fl9l with 



A 
A 2 



/ 4 
a (t) + -^u(r 



3ad 2 



A 4 



4r 3 

3ad 2 
" 4r 3 





ad 2 \ 




2r 3 J 




ad 2 \ 


(*+ 


2r 3 J 




^ 12 


-a'(r 
r 


) + — 



2a'fr 



ad" 



(21) 



(22) 



u"(r) + -u , (r)-4^(r)) I 



2u'(r 



9ad 2 u'(r) 6ad 



ad 2 



2r 4 



+ 



(23) 
(24) 
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The Dirichlet-Neumann problem for this part of the inhomogeneity is solved if the 
following differential equations hold (Note the additional minus sign since the right 
hand side is -^): 

4> B 2(r) = -A (a\cx",u) (25) 

B 2 '(r) + -B' 2 (r) - ±B 2 (r) + = -A 2 (a', at", u, u', u") (26) 

B» + lB4(r)-16^ = -A 4 (u,u',u") (27) 

We interpret the first two equations as defining equations for B2 and B4 given any a 
and u. The functions a and u are then determined by the third differential equation: 
We relate a and u in such a way that the third differential equation holds. In this case 
q (constructed from B2, B4) is a solution of the Dirichlet-Neumann problem. 
Carrying out the calculation explicitly we arrive at 

ad 2 + 2br 2 f41 , , 2br 2 - 7ad 2 m , , 55ad 2 - 34br 2 ... . 

24r ^ )(r) + 24r 2 ^ (t) + 24r^ ^ M 

85ad 2 -32br 2 „ , 245ad 2 , , 
U? ^--i2r^ (r) 

= - (ad 2 + 2br 2 ) a (6) (t) - ^ (-ad 2 + 1 Obr 2 ) a (5) (r) (28) 
ad2 ^ w (4) w J]5b_Mu(% 



3(8ad 2 -5br 2 ) f2 , 3ad 2 „ . 

What we have obtained is a differential equation relating the functions a and u. Can 
we find a gluing function a and a function u of bump-character related by d28i ? We 
interpret this equation as an inhomogeneous ordinary differential equation for u. The 
right hand side — independent of the choice of the gluing function — will be a function 
of bump-character. If we set initial conditions (j/ n '(1) = for n = 0,1,2,3 all 
derivatives of u will vanish at the point r = 1 by the properties of the right hand side. 
The only additional condition u has to satisfy is that the function and all its derivatives 
should also vanish at the point r = 2. Can we choose a right hand side, such that u and 
all its derivatives vanish also at r = 2? The answer is provided by the solution formula 
for inhomogeneous ordinary differential equations. We multiply (EHJ by aA i+2\, r i 
and obtain an equation of the following type 

u (4) (r) + f(r)u (3) (r) + gW^(r) + Mr)u'(r) + Mt)h(t) = JWM (29) 

Here f(r), g(r), h.(r),k(r) are well behaved, bounded, rational functions in the inter- 
val [1,2]. They also depend on the parameter b, which for notational simplicity we 
suppress to write out explicity most of the time. We will restore the argument b when 
it becomes important, as it will below. (The functions also depend on the other param- 
eters a, d. But this dependence does not interest us anyway.) We recall that 3~[<x] (r) is 
defined by the right hand side of d28b multiplied by ^i^bx 1 ¥ ® anc ^ nence depends 
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on the first to sixth derivatives of a. 3 All derivatives of 3~[a] (r) with respect to r vanish 
at the boundary r = 1 , 2. The dependence of CF on b is also notationally suppressed at 
this moment. 

We transform the system to a system of four ordinary differential equations of first 
order in the usual way: 



(30) 





/ 





1 





\ 






( ° 


\ 


u(r) 










1 





u(r) 


+ 







v(t) 













1 


v(t) 







\w(t)/ 


V 


Mr) 


-h(T) 


-gW 


-mj 


\w(r)J 









This we will write as 



XT 



A(r)x(r)|n(a[T)) . 



(31) 



with the obvious identifications. Derivatives are always taken with respect to r. We 
impose the boundary condition x(1) =0. 

The following standard arguments about ODEs are taken from Q. If <J>(r) is a 
fundamental matrix 4 of the homogeneous system, then 



XT 



0(t) 



<D _1 (t)n(a(t)) dt 



(32) 



is a solution of the inhomogeneous system satisfying the boundary condition x( 1 ) =0 
imposed above. It then automatically follows from d32l that x has altogether vanishing 
derivatives at r = 1 . We also need to impose the condition that x and all its derivatives 
vanishes at the other boundary at r = 2. How this can be achieved will be discussed 
next. But before we shall make the following observation: A priori our equations were 
real, so that O is a priori a real 4x4 matrix. However, it will turn out to be useful 
to complexify our system. Note that even if we use complex fundamental matrices 
we will get a real solution x of d32b . This follows from the fact that if O is some 
(complex) fundamental matrix, any other fundamental matrix is given by OC with 
some appropriate constant complex matrix C (compare Thm. 2.3 of 0). Obviously 
d32l is invariant under O i— ) (DC. 

Now we turn to imposing the second boundary condition: x(2) = 0. Again we 
observe that this boundary condition together with i32l implies that all derivatives, 
too, vanish at r = 2. Requiring two boundary conditions for an ODE of first order is 
clearly an overdetermination. However, in our case, the inhomogeneity in (13 It is not 
fixed. Hence we read x(r = 2) = as equation restricting the choice of the function <x. 
If we define v to be the vector whose components are the fourth column of the matrix 



, i.e. Vi 



<P, 4 1 , then the boundary condition is equivalent to 



r2 



v(t)J[<x](t)dt = 0. 



(33) 



-i 



Before we continue we remark that under a redefinition O i— > <DC we have (D 
C -1 ® -1 and hence v i— > C _1 v, showing that d33l is satisfied for one choice of O if 
and only if it is satisfied for any other. 

3 We exclude values b < — ad 2 /8. 

4 A matrix, whose n columns are the n linearly independent solutions of the homogeneous system. 
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To get an idea what may go wrong in trying to satisfy (I33t . we use the explicit 
expression for 3" and get, after some integrations by parts: 



i-2 



dr a'frl 



1 



r 2 v(b,r) + t 



ad 2 - 10br : 
ad 2 + 2br 2 



■v(b,r 



(4) 



4ad 2 - 5br 2 
ad 2 + 2br 2 

9 8ad 2 



v(b, r 



3 / 1 35br 2 - 28ad 2 
+ I - — — ^v(b,r 



r 2 ad 2 + 2br 2 



5br 2 

v(b,r 



r ad 2 + 2br 2 

72 ad 2 
r 5 " ad 2 + 2br 2 



(34) 



0. 



This clearly cannot be satisfied if one of the components functions in the square bracket 
is constant, since 



r2 



<x'(r)dr = a(2)-a(1) = 1. 



(35) 



More general, d34b cannot be satisfied if after left multiplication with some C _1 at least 
one of the component functions is constant. However, in the sequel we shall prove that 
at least for small parameters b this cannot happen. The idea is to prove this directly for 
b = and then use an inverse-function type argument to extend this to small b. Hence 
we first need to consider the case b = in some detail. 

3.6 The case b = 

In the case b = the system d28t reduces to 



Z (3) (r) + 55 // (r] 
r r z 



170 , 



490 



r J r 

lrV 6 '(r) + ra^{r) + 4ot^(r) - -ot^(r) + ^ 
2 r r 2 



ol [2] (t) 



72 



(36) 



a (/) (r) 



The homogeneous part is an equation of Euler-type, i.e. solved by simple powers r ni , 
so that the fundamental matrix is readily computed: 



0(r) 







r n 2 




r n 3 




r ri4 






Y 


(r n 2 


' 


( r n 3 


' 






(T n i; 


" 


(r n2 ) 


// 




// 




// 




III 


(r n 2) 


in 


(r n 3) 


in 


(r n 4) 


m) 



(37) 



The exponents n; e C are determined by the fourth-order polynomial obtained by 
inserting |x(r) = r nt into the left-hand side of d36l and equating this to zero. The 
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solutions are 



n-1 
n 2 



7 
2 
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n 3 = 2 + 



-72 + \/65]5) 3 
11 C 1 +iV3 



- + ( -72 + v 7 *^ 



1 -iV3) (-72 + ^6515 



n 4 



2 + 



2 (-72 + V6515 w ' 
11 C 1 - i y/3] 



1 +i\/3) (-72 + ^6515 



2 1-72 + ^6515 



(38) 
(39) 

(40) 

(41) 



Note that n-j,n.2 are real and 113,1x4 are complex conjugates. Of course we could 
construct two real solutions from the complex ones but, as already argued, we may 
continue to use a complex 4>. This is advisable since the complex form <37b . whose 
entries do not contain sums, has a simpler inverse. The fourth column of the inverse is 
the vector v. It reads: 



v r 





r 3 - 


n i 






(ui 


-rt 3 )(ni - 

r 3 - 


-rut 

n 2 


(m 


-n 2 ) 


(n 2 


-tv 3 )(tl 2 - 
r 3 - 


-TU 
n 3 


(n 2 


-ni) 


(n 3 


-TI4 ) (n 3 - 
r 3 - 


-ni 

n 4 


(n 3 
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(42) 



We can now explicitly evaluate the integral condition d3*3l . This leads to 



dr r 3 ~ ni f - lr 2 a (6) (r) + ra (5) (r) + 4a (4) (r) 



2 

28 



72 



— a (3) (r) + '-^cx {1 \t) 

T T A 



72 



(43) 







which is to be read as a condition for a, where the take the values computed above. 
Performing the partial integrations explicitly we can state the conditions in the form 



dr Wi(r)a'(r) =0 or 



r2 



-Wi(v 



(44) 



where the are four real, non constant, linearly independent functions on [1,2], 
whose exact form need not be spelled out here. Regarding our functions as elements 
of the Hilbert space L ( [1 , 2] , dr) this states that we need to choose a with prescribed 
orthogonal projections onto the four wt. This can always be achieved due to their 
linear independence (for linear dependent the projection conditions might turn out 
to be contradictory), while still leaving an infinite freedom for a. Hence we have 
shown that for b = the functions a and \x can be related such that the DN-problem 
is solved at first order. This we summarize as 
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Lemma: The DN-problem is solvable at the first order in the inverse glueing radius 
for the choice b = 0. In this case the glued and original data have the same overall 
ADM-energy. 



3.7 The case b > 

We now wish to extend the result to small positive b, which corresponds to a reduction 
of energy. The idea is to set up an implicit-function theorem. For this it is convenient to 
explicitly display the relevant dependencies on the parameter b, which we suppressed 
so far. 

Pick a glueing function ctoM which satisfies the four integral conditions (PHI . 
Define the Banach space of bump-functions 

S :={P € C°°([1,2]) | (3 (n) (1) = |3 (n) (2) =0 Vn e No} , (45) 

equipped with the max-norm on [1,2]. Next consider a map f : E x 25 — > E 4 , defined 
by: 



f : (b,(3) ' ) <D(b,r = 2) 



dtv(b,t)g-[b,a + P](t). (46) 



As already argued, this map indeed maps into E 4 even if we use complex fundamental 
matrices. Note also that the function ao(t) + |3(t) is a gluing function. The map d46t 
is a smooth map of Banach spaces, which is linear in oto + |3. 

We have shown above that f (0, 0) = 0. We wish to show that for any small b £ 
E there exists a (3 G 25 such that f (b, (3) = 0. This is ensured if the differential 
D 2 f (0, 0) : 25 -> E 4 is surjective. Since the map f is linear in oto + P, we arrive at 

r2 



drv(r) J[0,p](T), (47) 

i 



D 2 f(0,0) : P i — ) «D(0,t = 2) 

where v(0,r) = v(r) is explicitly given by the expression in (l42t . This indeed shows 
surjectivity and we have 

Theorem: For sufficiently small b € E the Dirichlet-Neumann-problem can be solved 
to first order in the inverse gluing radius. Hence the gluing can be used to reduce the 
ADM-energy. 



4 Conclusion 

Motivated by recent results concerning the radiation content at spacelike infinity of the 
most simple, time symmetric, conformally flat two-black-hole data (Brill-Lindquist 
data), we considered modifications of that data which render then Schwarzschild at 
infinity by using Corvino's gluing technique. We used Brill waves for the gluing and 
have given a perturbative argument (for large gluing radii) that this technique may be 
used to reduce the ADM-energy. We consider this as being a first small step towards 
making Corvino's technique more concrete, which also sheds some light onto the prob- 
lem of how to remove 'spurious' radiation at spacelike infinity. A natural extension 
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of this work would obviously consist in giving a rigorous proof for the possibility of 
energy reduction beyond finite-order approximation techniques. 

More generally, an interesting question to ask is what happens if one goes to 
smaller and smaller gluing radii. Could one reasonably expect this procedure to lead 
to a minimal overall ADM-energy for given hole-masses and -separations, so as to cor- 
respond to the ideal situation of no 'spurious' radiation? Or should one rather expect 
such a minimum not to exist? We do not know the answer to this question. 
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